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The spin-flip (SF) Eliashberg function is calculated from first-principles for ferromagnetic Ni to accu- 
rately establish the contribution of Elliott- Yafet electron-phonon SF scattering to Ni's femtosecond 
laser-driven demagnetization. This is used to compute the SF probability and demagnetization 
rate for laser-created thermalized as well as non-equilibrium electron distributions. Increased SF 
probabilities are found for thermalized electrons, but the induced demagnetization rate is extremely 
small. A larger demagnetization rate is obtained for non-equilibrium electron distributions, but its 
contribution is too small to account for femtosecond demagnetization. 
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Ultrafast demagnetization of ferromagnetic metals 
through excitation by a femtosecond laser pulse was dis- 
covered fifteen years ago by Beaurepaire et al. IJ. In 
spite of intensive investigations the microscopic origin of 
the ultrafast demagnetization could not be disclosed and 
continues to be controversially debated (see [2 for a re- 
cent review). Several mechanisms have been proposed to 
explain the observed ultrafast phenomenon [3l-fTT| . Most 
of these theories assume the existence of an ultrafast spin- 
flip (SF) channel, which would cause dissipation of spin 
angular momentum within a few hundred femtoseconds. 

Elliott- Yafet electron-phonon SF scattering has been 
proposed as a mechanism for ultrafast spin-dissipation 
[4j. Strong support in favor of electron-phonon medi- 
ated spin-flips as the actual mediator of the femtosec- 
ond demagnetization was made in a very recent work, 
in which ah initio calculated SF probabilities for ther- 
malized electrons compared favorably to SF probabili- 
ties derived from pump-probe demagnetization measure- 
ments |8j . While these results definitely favor the EUiott- 
Yafet SF scattering mechanism, the calculation of the 
electron-phonon scattering involved several serious ap- 
proximations. Applying the so-called Elliott approxima- 
tion |12j only spin-mixing due to spin-orbit coupling in 
the ah initio wavefunctions was included, but no electron- 
phonon matrix elements and no real phonon dispersion 
spectrum was considered. The thus-obtained SF proba- 
bility is however not a direct measure of demagnetization. 
Recent model simulations for thermalized hot electrons 
[9] using the Landau-Lifshitz-Bloch equation [13] and as- 
suming a fitted SF parameter did reproduce the exper- 
imental magnetization response, but couldn't assign the 
SF origin. Hence, it remains a crucial, open question 
whether laser-induced demagnetization can indeed be at- 
tributed to electron-phonon mediated SF scattering. 

Here we report an ah initio investigation to accu- 
rately establish the extent to which the Elliott- Yafet 
electron-phonon SF scattering contributes to fs demag- 



netization. To this end we perform ah initio calculations 
for ferromagnetic Ni, which ultrafast magnetization de- 
cay is well documented [J |51 [H]. We include the full 
electron-phonon matrix elements and phonon dispersions 
in our calculations. Introducing an energy-dependent 
SF Eliashberg function we compute SF probabilities and 
demagnetization rates for laser-heated thermalized elec- 
trons as well as laser-induced non-equilibrium electron 
distributions, from which we draw qualified conclusions 
on the possibility of phonon-mediated demagnetization. 

To treat phonon-mediated SF scattering at variable 
electron energies we define a generalized energy- and 
spin- dependent Eliashberg function, 

al,,F{E,^l) = ^Y. jj d^d]^gZi'n'{^) x 

5{UJ^. - M)5{Kn - mE^'n' - E) ,(1) 

which comprises initial and final electron states with 
quantum numbers kn, k'n' that interact through a 
phonon with frequency Q.=uic^u, v and q denote its mode 
and wavevector. M is the ionic mass, cr=t, i denote the 
spin majority, miniority components. For E—Ep (the 
Fermi energy) the SF part a'^^F{EF,fl) gives the SF 
Eliashberg function [TSl and the sum over all aa' corre- 
sponds to the standard Ehashberg function, a'^F{Ep, V,) 
[IB] . The (squared) electron-phonon matrix elements are 

gZxn'ici) = |Uq. • {nn\^RV\^U'n')\' , (2) 

where V is the potential, Uqi, the phonon polarization 
vector and are the eigenstatcs in the ferromagnet. 

Momentum conservation requires q = k' — k. SF scatter- 
ing becomes possible through the relativistic spin-orbit 
coupling. The majority, minority Bloch states |^'k„) and 
l^icri) ^® decomposed in pure spinor components 

l*L> = «L(o)+^L(")- l*L>-«L(i)+^L(o)' (3) 
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where the components are nonzero only if spin-orbit 
coupling is present and represent the degree of spin- 
mixing, which is a precondition for nonzero Skl^^k'n'- 

To study demagnetization we consider two quantities, 
SF probabilities and spin-resolved transition rates. The 
latter are defined as [11] 

^""'= J l'al^,F{E,n)UE){i-,u,{E+hn)) X 

{e{n)+N{n)) dndE. (4) 

Here N{^n) is the phononic Bose-Einstein distribution, 
the Fermi distribution, and Q(fl) the Heaviside function. 
Important for the effective demagnetization is the spin 
decreasing rate S~ , which corresponds to S^^, while the 
increasing one corresponds to 5"^^. 

An approximation of Eq. Q is helpful to achieve a 
faster evaluation and provide more insight in the pro- 
cess. Energy conservation during electron-phonon scat- 
tering requires i?k'n' — Eit^n — hfl, but the phonon energy 
hrt is usually very small (< 0.04 eV) compared to elec- 
tron related properties. Already in the standard Eliash- 
berg formulation Eq. ([I]) an energy difference between 
initial and final states is neglected while the ^-functions 
S{E^^ — E) are broadened with a parameter (0.03 eV, 
here). Similarly, one can neglect the energy variation 
due to hn, in the Fermi function fcr{E + hfl), as long as 
the temperature is high enough. Wc can then rewrite 
spin-resolved transition rates in the form 

5""' = I w^,,{E)UE){l - U'{E))dE, (5) 

where we introduced the energy- and spin- dependent 
specific scattering rate for electrons Wo-o-' given by 

Waa'iE)^ / dnal,,F{E,n){l+2N{n)). (6) 
Jo 

Note that w^i{E) — wi^{E). All calculations were 
checked against a more accurate numeric implementation 
not involving this approximation. The SF probability for 
an electron with energy E is defined as the ratio of the 
SF part to the corresponding total counterpart, ps{E) = 
2w-^l{E) I Ylicra' "^crcr' {E) . Aualogously, the total SF prob- 
ability during a scattering event can be defined as 

Ps^{,S- + S^)lY^S--' . (7) 

Although the SF probability has been used in recent dis- 
cussions of laser induced-demagnetization [51 [TS], it is 
actually not the crucial quantity (as a high but equal SF 
probability for both spin channels would not cause a de- 
magnetization) . We define therefore the normalized de- 
magnetization ratio, Dg = {S^ ^ I Tlaa' ^'^'^ ' which 
tracks the difference of magnetic moment increasing and 
decreasing SF contributions. 



To investigate phonon-induced demagnetization in 
laser-excited Ni we proceed now in three steps. First, we 
compute the ab initio SF probability Ps for equilibrium 
Ni, i.e., for E = Ep. Second, we compute SF probabihties 
Ps for laser-heated Ni, by treating a range of electron en- 
ergies that correspond to those in a hot, thermalized elec- 
tron gas after laser-excitation. Thermalization to elec- 
tron temperatures Tg of a few thousand K occurs quickly 
within about 200 fs after the laser pulse, but the hot elec- 
trons are not in equilibrium with the lattice and the lat- 
tice temperature is not altered significantly. In the third 
step we consider the SF probability for non-equlibrium 
(NEQ) electron distributions [19j that are expected to be 
present within ~100 fs after laser stimulation. Demag- 
netization ratios Dg are subsequently evaluated for these 
three situations. The results obtained in these steps are 
furthermore compared to values which we compute with 
the so-called Elliott relation (see below). 

An ah initio evaluation of the SF probability of equi- 
librium Ni requires calculated phonon dispersions and 
a relativistic electronic structure. Such calculation has 
previously been done for paramagnetic Al |15j . but has 
not yet been accomplished for ferromagnets. An approx- 
imation was introduced years ago by Elliott 12J, who 
pointed out a possible source of SF scattering arising 
from the spin-mixing of eigenstates. Employing several 
assumptions, viz. a paramagnetic metal, nearly constant 
electron-phonon matrix elements, 6kn constant in the 
Brillouin zone, and ^ a^„, Elliot derived a rela- 
tion between the spin lifetime rg for a general kind of 
scattering event with lifetime r. This so-called Elliott 
relation uses the Fermi surface averaged spin-mixing of 
eigenstates (fe^) ^ Y.a.n S dV\Kn?SiKn " Ef) and pre- 
dicts the SF probability Pf = [ts/t)-^ = 4(52). 

In a similar way as introduced above, the influence of 
spin-mixing on the SF probability in laser-heated Ni can 
be evaluated. We define a SF density of states (DOS) as 

nn{E)^Y. h^\^^n?5{El^-E)- (8) 

A generalized Elliott SF probability for an electron with 
energy E is then given as Pg (E) = 4n^^(i?)/ n{E) (with 
n{E) the total DOS) which yields the standard Elliott ex- 
pression (6^) in the limit <C and E = Ef- The to- 
tal SF probability Pg of a laser-heated system with elec- 
tron distribution fa{E) is obtained from Eqs. ([t]) and ([5|, 
where w^i{E) is replaced by n^^{E) and w{E) by n{E). 
Note that although the treatment is intended for phonon 
scattering the Elliott relation in fact does not take the 
character of scattering involved into account. Also, the 
assumption of a paramagnetic material is essential in El- 
liott's derivation as this permits SF scattering in each k 
point in the spin-degenerate majority, minority bands at 
Ep- Experimentally the Elliott relation was found to be 
valid up to a multiplication by a material specific con- 
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FIG. 1. (Color online) Ab initio calculated Eliashberg 
a'^F{EF,0.) and SF Eliashberg a'^^^FiEp,^) functions of Ni 
in equilibrium. 

stant with variation smaller than one order of magnitude 
for various paramagnetic metals i20j . Recently it has also 
been applied to ferromagnetic metals O [18] , even though 
for exchange-split ferromagnetic bands there exist far less 
k points at which spin-degenerate bandcrossings occur. 

We have tested the implementation by computing first 
Al and Ni in equilibrium at low temperature (<300K). 
Our calculations are based on the density functional the- 
ory (DFT) within the local spin-density approximation 
(LSDA), see [H] for details. For Al our calculated a'^^F 
is of the order of 10^ smaller than a^F and in good agree- 
ment with the existing previous result [TS] . The ab initio 
calculated SF and non-SF Eliashberg functions of equi- 
librium Ni are shown in Fig. [l] For Ni the computed SF 
a^^F function is only about 50 times smaller than the or- 
dinary a^F function; this is due to the larger spin-orbit 
coupling. The resulting total SF probability, Ps=0.04, is 
given in Table |l] To estimate the accuracy of the Elliott 
approximation we have calculated the Elliott SF proba- 
bility and obtain Pg =0.07. This value is in rough agree- 
ment with P| =0.10 computed in Ref. [H]. Thus wc find 
that the Elliott relation overestimates the SF probability 
in equilibrium Ni by about a factor two. 



TABLE I. Calculated spin-flip probabilities Ps and demagne- 
tization ratios Ds for laser-pumped Ni. Results are given for 
equilibrium (low T), for thermalized electrons at a high Fermi 
temperature Te, and for the non-equilibrium (NEQ) electron 
distribution created by fs laser-excitation. Results obtained 
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for the approximate Elliott SF probability Pg (this work and 
|18p are given for comparison. 
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FIG. 2. (Color online) Energy-resolved electron-phonon total 
and SF scattering rates w^E) and w-^ilE) of Ni, and normal- 
ized SF probability Ps{E) and approximate SF probability 

,2 

Ps (E) obtained from the Elliott relation. 



Next we turn to the topic of current controversy, 
the actual amount of phonon-induced demagnetization 
in laser-excited Ni. In Fig. [2][top) we show calculated 
energy- resolved SF and non-SF scattering rates {w^i{E) 
and w{E)). Note the strong energy variations of w{E). 
In Fig. [2j^bottom) we compare the computed electron- 
phonon SF probability Ps{E) to that obtained from the 
Elliott relation. At some energies, e.g., 0.5-1 eV, these 
two quantities are nearly the same, but at other energies 
there is no direct relation other than that SF probability 
is large where band states are present. An interesting dif- 
ference in the context of ultrafast demagnetization is the 
suppression of Ps{E) around Ep, which is not captured 
by Pg (E). The features of Ps{E) that are not captured 
by Pf (E) can be understood by comparing Eqs. ^ and 
([S]). One of the differences is the presence/absence of 
summation over destination eigenstates k'n'. The latter 
are restricted in Eq. tin by the construction of ff^J^^^/^/ 
to correspond to a different spin than the source state 
kn. The number of available end states is however not 
taken into account in Elliott formula (which, derived for 
a paramagnetic metal, assumes that the same number 
of states is available for both spins, and hence suppresses 
this distinction). The mentioned discrepancy between Ps 
and Pg above Ep is thus easily explained by the lack of 
states with the same energy and opposite spin in the Ni 
DOS (see Fig. |3|. Hence, the Elliott relation fails for 
ferromagnets in strongly exchange-split energy regions. 

After laser-excitation electrons equilibrate quickly due 
to electron-electron scattering at a high electron tem- 
perature Te of the order of thousands K. To describe 
this situation we use appropriate fa-{E), but note that 
the chemical potential must be adjusted also. Spin con- 
servation leads to differences between f^{E) and fi{E), 
namely fi{E) has a lower chemical potential than f-\{E) 
in Ni due to the shape of its DOS. SF probabilities Ps 
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FIG. 3. (Color online) Spin-resolved DOS (filled areas) and 
phonon induced spin-flips (arrows) of NEQ and electron ther- 
malized Ni. The equilibrium DOS is shown by thin lines. SF 
transitions are significantly different at energies above and 
below Ef (=0 eV). The arrows thickness corresponds to the 
transition rate, its direction and length give which direction is 
dominant and how much. The amount of laser redistributed 
electrons has been enlarged to improve visibility. 



computed for several are given in Table |T] With in- 
creasing Te Pg increases, too. Also the Elliott SF proba- 
bility Pg increases with Te, but it deviates still from P5. 
A previous work [B] used a Gaussian smearing to stimu- 
late a thermalized system (vi^ithout Ep adjustment) and 
obtained Pg saO.18. Our values are smaller, but note 
that the way the thermalized distribution is described is 
different. 

As mentioned before, a large SF probability does not 
necessarily imply a large demagnetization. Evaluating 
the demagnetization rate dM/dt=2iiB{S^—S^) for ther- 
malized electron distributions we obtain quite small val- 
ues, of the order of O.OS/is/ps. The reason is that not 
just a large SF probability, but also an imbalance between 
f\{E) and fi{E) is essential for a magnetization change. 
The distributions of spin populations specific to Ni imply 
that for thermalized electrons below Ep most spin-flips 
increase the spin moment, spin-reducing transitions oc- 
cur only above Ep. In that region the SF scattering rate 
is however very low (Fig.[2|. The situation is illustrated 
in Fig. [3j As a consequence the spin-decreasing rate 
(S'"-S'+) is thus much lower than the SF rate {S~+S+), 
and in addition it exhibits only a weak temperature de- 
pendence. Hence we find that phonon-mediated SF scat- 
tering in thermalized Ni cannot be the mechanism of the 
observed ultrafast demagnetization. 

One remaining possibility for a fast demagnetization 
is an enhanced SF rate in the NEQ distribution present 
immediately after the laser pulse. Previous ab initio cal- 
culations showed that minority-spin electrons are excited 
more than majority-spin ones, see |19) . Assuming a 1.5- 
eV pump-laser and a simplified step-like electron distri- 
bution reduced by about 5% in the 1.5-eV energy win- 
dow below Ep, the calculated demagnetization ratio Ds 
is higher than for thermalized distributions (Table |l]). A 



critical role is played here by holes deep below Ep with 
high SF probability as well as a significant difference be- 
tween majority and minority occupations (see Fig.js]). An 
important yet unknown element in estimating the demag- 
netization is the laser fluence. Nonetheless, we find that 
phonon-mediated demagnetization in Ni is much more 
effective in the NEQ state than in the thermalized state, 
as was proposed recently for Gd An important as- 
pect is the time scale on which the NEQ demagnetiza- 
tion is active. Electron thermalization proceeds fast in 
Ni and transforms the initial NEQ distribution to a ther- 
malized one in ^200 fs. A rough estimate of the de- 
magnetization in this time-window is 0.1/iB, i-e. smaller 
than the observed experimental demagnetization. The 
precise amount of the demagnetization depends however 
on the time-evolution of the distributions, which requires 
further investigations. 

Using relativistic ab initio calculations we have evalu- 
ated the phonon-induced SF probability and demagneti- 
zation in laser-pumped Ni. A strong dependence of these 
quantities on the electron energy is observed, which is 
not tracked by the Elliott approximation. In the electron 
thermalized state Elliott- Yafet phonon-mediated demag- 
netization is too small to explain the ultrafast demagneti- 
zation, despite reasonably large SF probabilities. We find 
that Elliott- Yafet SF scattering contributes more to the 
demagnetization for NEQ distributions immediately after 
the fs laser-excitation. We note lastly that the existence 
of other fast SF channels [SHZl |TT] cannot be excluded. 
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